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>— P ' Abstract. Near-infrared spectroscopy is used to study the kinematics and excitation mechanisms of H2 and [Fe 11] 

gas in a sample of mostly Seyfert 1 galaxies. The spectral coverage allows simultaneous observation of the JHK 
bands, thus eliminating the aperture and seeing effects that have usually plagued previous works. The H2 lines 
are unresolved in all objects in which they were detected while the [Fe II] lines have widths implying gas velocities 
^— V , of up to 650 km s~^. This suggests that, very likely, the II2 and [Fell] emission does not originate from the same 

■~~^ ■ parcel of gas. Molecular II2 lines were detected in 90% of the sample, including PG objects, indicating detectavel 

1-^ ' amounts of molecular material even in objects with low levels of circumnuclear starburst activity. Analysis of the 

observations favors thermal excitation mechanisms for the II2 lines. Indeed, in NGC3227, Mrk766, NGC4051 and 
O ' NGC4151, the molecular emission is found to be purely thermal but with heating processes which vary among 

i-i ' the objects. Thermal excitation is also confirmed by the rather similar vibrational and rotational temperatures 

C/3 , in the objects for which they were possible to derive. [Fen] lines are detected in all of the sample AGN. The 

Cu ' [Fen] 1.254/im/Pa/3 ratio is compatible with excitation of the [Fen] lines by the active nucleus in most Seyfert 1 

galaxies, but in Mrk766 the ratio implies a stellar origin. A correlation between H2/Br7 and [Fen]/Pa/3 is found 
for our sample objects supplemented with data from the literature. The correlation of these line ratios is a useful 
diagnostic tool in the NIR to separate emitting line objects by their level of nuclear activity. X-ray excitation 
models are able to explain the observed H2 and part of the [Fe n] emission but fails to explain the observations in 
Seyfert 2 galaxies. Most likely, a combination of X-ray heating, shocks driven by the radio jet, and circumnuclear 
star formation contributes, in different proportions, to the H2 and [Fen] lines observed. In most of our sample 
objects, the [Fen] 1.257/im/1.644/im ratio is found to be 30% lower than the intrinsic value based on current 
atomic data. This implies either than the extinction towards the [Fen] emitting clouds is very similar in most 
objects or there are possible inaccuracies in the A-values in the Fen transitions. 

Key words, galaxies: Seyfert - infrared: spectra - molecular processes infrared: galaxies - line: formation 

1. Introduction or shocks from supernova remnants). This ambiguity is 

most evident when interpreting the spectral lines of low- 
One of the fundamental problems in active galactic nu- ionization species such as [Fen] and molecular hydro- 
clei (AGN) research is to determine the dominant exci- g^^^^ ^^^^ g^^g of ji^^eg ^^^ detected in galaxies display- 
tation mechanisms of the narrow line emitting gas, i.e. ^^^ varying degrees of nuclear activity, from objects clas- 
whether it is due to non-stellar processes (e.g. photoion- gjggj as starburst-dominated to those classified as AGN- 
ization from a central source or shocks from a radio jet) dominated. Do these lines are produced by a unique mech- 
or to stellar processes (e.g. photoionization from OB stars ^nism in the objects where they are observed, evidenc- 
ing the so-called AGN-starburst connection? or are they 



Send offprint requests to: A. Rodn'guez-Ardila 



Alberto Rodriguez- Ardila et al.: H2 and [Fen] in AGN 



formed by different mechanisms (or some combination 
thereof) in different classes of objects? 

Transitions within the ground state of H2 are often 
bright in near-infrared observations of AGN. The origin 
of these lines has been highly debated, with no definitive 
consensus yet reached. Three types of excitation mech- 
anisms may produce the H2 lines: (a) UV fluorescence 
(e.g Black & Van Dishoek lTMTjl : (b) shocks (Hollenbach 
& McKee [l989|l : and (c) X-ray illumination (Maloney et 
al. I1996|l . UV fluorescence is usually regarded as a non- 
thermal excitation process and shocks and X-rays as ther- 
mal ones. Theoretically, each of these processes produces 
a different spectrum, and the relative intensities between 
emission lines of H2 can be used to discriminate the dom- 
inant excitation mechanism. So far, however, the results 
are controversial. For example, Veilleux et al. H1997|l . claim 
that shocks associated with nuclear outflows are the likely 
source of H2 in AGN while Quillen et al. H1999|l relate 
the H2 emission to local star formation. Very recently, 
Rigopoulou et al. (j2()()2(l carried out a survey of molec- 
ular hydrogen emission using a sample of starburst and 
Seyfert galaxies observed by ISO. They conclude that the 
H2 emission in Seyferts is most likely a combination of that 
from photodissociation regions, shock-heated gas, and X- 
ray heated gas, thus favoring thermal excitation over the 
non- thermal. However, most of the above conclusions have 
been obtained using a small number of objects with a clear 
bias towards Seyfert 2. In addition, the aperture used in 
the ISO observations is usually large which allows the con- 
tamination of the nuclear spectra by the host galaxy. 

Similar problems, namely production mechanisms and 
location of the emitting region, are associated with the 
interpretation of [Fe 11] emission in AGN. Forbes & Ward 
H1993|l found a tight correlation between [Fell] and the 
6-cm radio emission for both starburst (SB) and Seyfert 
galaxies. This was taken as strong evidence of shock exci- 
tation from the radio jets and/or SNRs for the production 
of [Fell] in AGN. However, Simpson et al. I|1996|l argued 
that photoionization by the central engine must be the 
dominant excitation mechanism. More recently, Mouri et 
al. H2()()()|l support this scenario by proposing that heating 
by X-rays is more important in the production of [Fe ll] 
and that photoionization by a central source or by shocks 
can be distinguished by the electron temperature of the 
[Fe 11] region: 8000 K in the former case and 6000 K for the 
later. Another useful indicator proposed to distinguish be- 
tween SB and Seyfert activity for the origin of [Fe 11] is the 
ratio [Fell] 1.257^m/Pa/3 (Goodrich et al. 1994, Alonso- 
Herrero et al. 19971 . Observationally, there is an increas- 
ing progression in that ratio from pure photoionization to 
pure shock excitation. In SB galaxies it is observed to be 
in the range of 0.1-0.5 (Simpson et al. I1996|l . in reason- 
able agreement with the prediction of SB models of Colina 
P993|l . while in AGN it is expected to be in the range of 
0.6-3 if heating by X-rays dominates. 

In order to address the above problems and to increase 
our knowledge about the origin of the H2 and [Fen] lines, 
we have observed a sample of 19 Seyfert 1, 2 Seyfert 2 and 



1 starburst galaxies. These objects are part of an ongoing 
program aimed at studying the near-infrared properties 
of AGN. For most of our target objects, we provide the 
first reports of molecular and [Fe 11] emission. The remain- 
ing paper is organized in the following manner: §2 presents 
the observations and data reduction; §3 discusses the kine- 
matics of the molecular and [Fe 11] gas inferred from their 
line profiles; §4 discusses the origin of the H2 lines; §5 deals 
with the emission mechanisms leading to the [Fe 11] lines; 
concluding remarks are given in §6. 



2. Observations and data reduction 

Near-infrared (NIR) spectra from 0.8-2.4 /im were ob- 
tained at the NASA 3 m Infrared Telescope Facility 
(IRTF) on April 20-25/2002 with the SpeX spectrom- 
eter (Rayner et al. 2003). The detector consisted of a 
1024x1024 ALADDIN 3 InSb array with a spatial scale 
of 0.15" /pixel. Simultaneous wavelength coverage was ob- 
tained by means of prism cross-dispersers. A 0.8" x 15" slit 
was used during the observations, giving a spectral reso- 
lution of 360 km s~^. For the objects in which the host 
galaxy is clearly detected on the DSS images, the slit was 
oriented along the major axis of the target. Otherwise, 
it was aligned to the parallactic angle in order to min- 
imize slit loses because of differential refraction. During 
the different nights, the seeing varied between 0.7"-l". 
Observations were done nodding in an ABBA source pat- 
tern with typical individual integration times of 120 s and 
total on-source integration times between 30 and 50 min- 
utes. For some sources, we took spectra on subsequent 
nights which were, after reduction, combined to form a 
single spectrum. During the observations, AO V stars were 
observed near each target to provide telluric standards at 
similar air masses. They were also used to fiux calibrate 
the sample. Table^shows the log of the observations and 
the atmospheric conditions during each night. The galax- 
ies are listed in order of right ascension. 

The spectral extraction and wavelength calibration 
procedures were performed using SPEXTOOL, the in- 
house software developed and provided by the SpeX team 
for the IRTF community (Gushing et al. 120031 Vacca et 
al. 12003(1 ^. Each pair of AB observations was reduced 
individually and the results averaged to provide a final 
spectrum. A 1.6"— 2" aperture, depending on the seeing 
conditions during the observations, was used to extract the 
spectra. Observations of an argon arc lamp enabled wave- 
length calibration of the data; the rms of the dispersion 
solution was, on average, 0.18 A. No effort was made to 
extract spectra at positions different from the nuclear re- 
gion even though some objects show evidence of extended 
emission. The small beam employed (usually 1.8" x 0.8") 
and the good seeing during the observations makes this 
data set truly representative of the nuclear emission, with 



^ SPEXTOOL is available from the IRTF web site at 
_http://irtf. ifa.hawaii.edu/Facility/spex/spex.htnil_ 
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Table 1. Observation log and basic galactic properties for the sample. The objects are listed in order of right ascension. 
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Integration 




PA 


r" 


ID 


Galaxy 


Type 


z 


E(B-V)g 


Observation 


time (s) 


Airmass 


D 


(pc) 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


(8) 


(9) 


(10) 


1 


Mrk 1210 


Sy2 


0.01406 


0.030 


2002 Apr 25 


2700 


1.25 


58 


220 


2 


Mrk 124 


NLSl 


0.05710 


0.015 


2002 Apr 23 


2640 


1.16 


10 


990 


3 


Mrk 1239 


NLSl 


0.01927 


0.065 


2002 Apr 21 
2002 Apr 23 


1920 
1920 


1.08 
1.15 


0.0 
0.0 


335 


4 


NGC3227 


Syl 


0.00386 


0.023 


2002 Apr 21 
2002 Apr 25 


720 
1080 


1.00 
1.02 


158 
158 


67 


5 


H1143-192 


Syl 


0.03330 


0.039 


2002 Apr 21 


1920 


1.31 


45 


520 


6 


NGC3310 


SB 


0.00357 


0.022 


2002 Apr 21 


840 


1.21 


158 


56 


7 


NGC4051 


NLSl 


0.00234 


0.013 


2002 Apr 20 


1560 


1.17 


132 


37 


8 


NGC4151 


Syl 


0.00345 


0.028 


2002 Apr 23 


1800 


1.10 


130 


58 


9 


Mrk 766 


NLSl 


0.01330 


0.020 


2002 Apr 21 
2002 Apr 25 


1680 
1080 


1.06 
1.02 


112 
112 


230 


10 


NGC4748 


NLSl 


0.01417 


0.052 


2002 Apr 21 
2002 Apr 25 


1680 
1440 


1.29 
1.21 


36 
36 


254 


11 


Mrk 279 


NLSl 


0.03068 


0.016 


2002 Apr 24 


3600 


1.54 


0.0 


480 


12 


NGC 5548 


Syl 


0.01717 


0.020 


2002 Apr 23 


1920 


1.05 


112 


298 


13 


Mrk 478 


NLSl 


0.07760 


0.014 


2002 Apr 20 


3240 


1.06 


0.0 


1200 


14 


NGC 5728 


Sy2 


0.01003 


0.101 


2002 Apr 21 


960 


1.31 


36 


160 


15 


PG 1448+273 


QSO 


0.06522 


0.029 


2002 Apr 24 


2160 


1.01 


108 


1020 


16 


Mrk 291 


NLSl 


0.03519 


0.038 


2002 Apr 21 


2520 


1.04 


84 


550 


17 


Mrk 493 


NLSl 


0.03183 


0.025 


2002 Apr 20 
2002 Apr 25 


1800 
900 


1.07 
1.04 


0.0 
0.0 


500 


18 


PG 1612+261 


QSO 


0.13096 


0.054 


2002 Apr 23 


2520 


1.10 


107 


2050 


19 


Mrk 504 


NLSl 


0.03629 


0.050 


2002 Apr 21 


2100 


1.04 


138 


570 


20 


IH 1934-063 


NLSl 


0.01059 


0.293 


2002 Apr 21 
2002 Apr 25 


1200 
720 


1.13 
1.17 


150 
150 


170 


21 


Mrk 896 


NLSl 


0.02678 


0.045 


2002 Apr 23 
2002 Apr 24 
2002 Apr 25 


1440 
1200 
1200 


1.21 
1.18 
1.17 


150 
150 
150 


420 


22 


Ark 564 


NLSl 


0.02468 


0.060 


2000 Oct 10 


1500 


1.05 


0.0 


390 



Radius of the integrated region. 



little contamination of the host galaxy or the extended 
narrow line region. 

Final wavelength-calibrated target spectra were di- 
vided by the spectrum of an AOV star observed at a similar 
airmass to eliminate telluric contamination. Typical air- 
mass differences between the target and the standard were 
0.05-0.1. However, despite these small differences, a small 
residual, mainly in the interval 1.995/im— 2.087/im, was 
found in some objects due to bad telluric band cancella- 
tion. After division by the AOV standard, each spectrum 
was multiplied by a blackbody of temperature equal to 
the effective temperature of the star to restore the true 
continuum shape of the target. 

Flux calibration was carried out by normalizing to the 
K-band magnitude of the corresponding telluric standard 
spectrum. We estimate that the uncertainty in this pro- 
cedure is ^10% by comparing our calibration with the 
2MASS fluxes available in the literature. The agreement 
in the continuum flux level in the overlap region for the 
different dispersion orders was excellent with errors less 
than 1%. The spectra were then corrected for redshift, de- 



termined from the average < z > measured from the po- 
sitions of [Siii] 0.953^m, PaJ, Hei 1.083^m, Oi 1.128/im, 
Pa/3 and Br7. Finally, a Galactic extinction correction 
as determined from the COBE/IRAS infrared maps of 
Schlegel, Finkbeiner & Davis (1998) was applied. The 
value of the Galactic E(B- V) used for each galaxy is listed 
in Col. 5 of TableQ] Final reduced spectra, in the spectral 
regions of interest in this work, are plotted in Figures ^ 
toEl 



3. Kinematics of the H2 and [Feii] lines 

The large spectral coverage and medium spectral resolu- 
tion of our data (360 km s^^) allow us to discuss how the 
widths of the H2 lines compares to that of other narrow 
lines. The goal of this analysis is to check if the molecu- 
lar hydrogen is kinematically linked to the NLR gas, set- 
ting important constrains on its location and origin. To 
our knowledge, no previous results about this issue are 
available in the literature, mainly because of the lack of 
adequate spectral resolution in the NIR. Recent NICMOS 
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24000 



Fig. 1. Final reduced spectra, in the Earth's frame, centred in near Pa/3 (left panel), the H band (1.60/L/,ni, middle 
panel), and Br7 (right panel). The abscissa is the observed flux in units of 10"'^^ erg cm^^ s~^ A~^. The identified 
[Fell] (dashed lines) and H2 lines (dotted lines) are marked in the spectra. Note that in the H band, no H2 lines 
are detected except in NGC5728 and NGC3227 where (1-0)S(7) 1.747/zm is visible. We show, however, the expected 
position of two of the strongest lines predicted by the Black & van Dishoeck's (1987) models. 



imaging of H2 emission by Quillen et al. I|1999|l for a sam- 
ple of Seyfert galaxies found H2 on scales of a few hundred 
parsecs from the nucleus for three objects (NGC 5643, 
NGC2110 and Mrkl066). The H2 emission was found to 
be coincident with [Olli] and Ha-|-[Nii], suggesting that 
the molecular gas may follow the NLR gas. We can test 
this connection by examining the the line profiles: if the 
line widths of the molecular gas are similar to that of 
other narrow emission lines, then, under the assumption 
that the widths of the NLR lines reflect the large scale 
motions of the emitting clouds in the gravitational poten- 
tial of the central mass, similar widths of molecular and 
forbidden lines suggest that they are co-spatial. 

TableEllists the FWHM in km s^^ measured in [Siii] 
0.9531 ^nJi, [Sivi] 1.963 /im, [Fen] 1.2567 ^i^Ji and H2 
(1,0)S(1) 2.121 fim. These values were obtained by fitting 
a Gaussian to the observed profile and then adopting the 



FHWM of the fit as representative of the FWHM of the 
line. The software LINER (Pogge & Owen HMSIl was used 
for this purpose. Afterwards, we corrected the measure- 
ments for instrumental resolution. In all cases, the lines 
in Table |2] were well represented by a single Gaussian 
component. For [Siii], [Fen] and H2 (1,0)S(1) 2.121/xm, 
no statistically significant shifts relative to the fiducial z 
were found while the centroid of [Si Vi] was observed to be 
blueshifted by up to a few hundreds km s~^ in some ob- 
jects, as is typical for coronal lines (Penston et al. I1984|l . 



Because this latter line is out of the main focus of this pa- 
per, we will not make further comments about its shifts. 
Typical errors, based on the uncertainties in the place- 
ment of the continuum, are ~30 km s~^. We have chosen 
II2 (1,0)S(1) 2.121/im as representative of molecular hy- 
drogen because it is one of the strongest H2 line in the 
NIR. In addition, it is isolated from other emission fea- 
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Fig. 2. Same as Fig. ^ 



tures, minimizing the effects of blending. Note that the the 
width and shape of other H2 lines within a given galaxy 
is rather similar to that of H2 2.121/im. No Hi lines were 
used in this analysis because, for more than half of the 
sample, the deblending of the observed profile into the 
components emitted by the NLR and BLR is highly un- 
certain. 

The results show that in all objects the width of H2 
(1,0)S(1) is at the limit of the spectral resolution (360 
km s~^) while [Fen] is resolved in 50% of the sam- 
ple. [Siii] 0.953 urn (IP==23.3 eV), considered the NIK 
counterpart of [Oiii] 5007 A, is unresolved in 60% of 
the sample, and in the remaining 40%, the trend is to 
have similar or larger width than [Fen]. The coronal line 
[Sivi] 1.963/im (IP=167 eV), when resolved, is signifi- 
cantly broader than H2 (1,0)S(1) and, except for a few 
cases, is also broader than [Fen]. Overall, in those objects 
in which the forbidden lines were resolved, we found in- 
creasing line widths with increasing ionization potential 
of the emitting species. A similar trend using optical and 
NIR lines was found by Rodriguez- Ardila et al. ( 2002|) . 



They showed, for a small sample of Seyfert 1 galaxies, a 
positive correlation between FWHM and ionization poten- 
tial. In the optical region, this type of correlation has been 
firmly established for a large number of AGN (Penston et 
al. HHHl De Robertis & Osterbrock [TUM Evans [T^55jl . 
This trend is usually interpreted as a stratification of 
the emitting regions in the sense that the highest ion- 
ized gas is located closer to the central source than low- 
ionization/neutral NLR gas. 

NGC3227, NGC4151 and PG 16124-261 are counter- 
examples to the above trend. In these objects, [Fen] is 
the broadest NLR line, even when compared to high ion- 
ization lines such as [Six] 1.252/im (IP=328 eV). From a 
kinematical point of view, the bulk of the [Fe n] and II2 
emission cannot originate in the same volume of gas. Part 
of [Fen] must arise from an additional source. We sug- 
gest that this extra source may be associated with shock 
excitation from the radio jet. In fact, Knop et al. H1996|l 
proposed this scenario for NGC4151. They found [Fen] 
broader than the narrow component of Pa/3 and [S ix] , 
the latter two lines also unresolved in their nuclear spec- 
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Fig. 3. Same as Fig. ^ 



trum. Knop's FWHM for [Fell] in NGC4151 is close to 
ours (434 km s^^ vs 505 km s^^, respectively). 

The smaller values of the FWHM of H2 relative to 
other narrow lines may suggest that the molecular gas is 
concentrated in the external/extended NLR or possibly 
even in the host galaxy, where the gravitational effects of 
the black hole on the emitting gas are lower than in the 
BLR. This scenario seems unlikely, however, because for 
the large majority of the objects the size of the integrated 
region covered by the spectra is no larger than 500 pc (see 
Col. 10 of Table ^. The alternative is to consider that the 
molecular gas is not kinematically linked to the NLR gas 
even though the two are co-spatial. 

NGC 3227 supports this scenario. Figure ^ji shows 
a comparison of the line profiles of H2 2.121^m and 
[Fen] 1.257/im. It can be seen that the latter line is much 
broader than the former although both are emitted from 
a region ~70 pc in radius. Note that [Fell] 1.257/im is 
slightly asymmetric towards the blue, mainly in the wing 
of the line, but it is not possible to tell from the data if the 
asymmetry is related to [Fe 11] or it is due to He i 1.252/xm. 



NICMOS imaging in H2 (l-O)S(l) 2.121/im carried out by 
Quillen et al. ( 1999) on this object shows that the molecu- 
lar hydrogen emission originates in a 100 pc diameter disk 
that is not directly associated with either the [O ill] emis- 
sion or the 18 cm radio emission. Moreover, Rodriguez- 
Ardila & Viegas §20031) detected 3.3/um PAH emission in 
the inner 120 pc of this object, implying circumnuclear 
star formation within that region probably related to the 
disk of molecular gas. The differences in line widths re- 
ported here supports the view that both the NLR and the 
molecular gas are co-spatial but emitted by different vol- 
umes of gas and most probably originating from different 
excitation mechanisms. 

Mrkl210 is another interesting case for a kinematical 
decoupling of the molecular gas. The forbidden NLR lines 
are remarkably similar in FWHM (see Table EJ; however, 
H2 (l-O)S(l) 2.121/im is significantly narrower, almost half 
the FWHM value of the [Fell] 1.257/im (FigureEti). Radio 
observations in the CO(l^O) 115 GHz and C0(2^1) 230 
GHz carried out on this object (Raluy et al. ll998|l show a 
high level of star forming activity within the inner 5 Kpc. 
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Table 2. FWHM (in km s^^), corrected for instrumental resolution, for H2 and other important lines measured in 
spectra sample. 





[Siii] 


[Si VI] 


[Fen] 


H2 (1,0)S(1) 


Oi 


Galaxy 


0.9530Mm 


1.963/im 


1.257Mm 


2.121^m 


1.128/im 


Mrk 1210 


685 


660 


650 


380 




Mrk 124 


420 


360 


420 


395 


1330 


Mrk 1239 


865 




456 




1160 


NGC3227 


470 


420 


560 


360 


3110 


H 1143-192 


360 


556 


404 




1630 


NGC3310 


360 




360 
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NGC4051 


360 
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3520 
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360 


590 


440 
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1550 
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440 


370 


2790 


NGC 5548 


360 


480 
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Fig. 4. Comparison of the line profiles of [Fell] 1.257/im 
(full histogram) and H2 2.122/im (dashed histogram) for 
(a) NGC 3227 and (b) Mrkl210, for which the observed 
profile of Br7 (solid line) has also been plotted. 



Moreover, this Seyfert 2 has a strong H2O megamaser of 
nuclear origin, arising from the innermost parsecs of the 

mm 



host galaxy (Braatz et al. 11994(1 . Because our integrated 
spectrum covers a region of '^200 pc in radius, it is likely 
that it includes a significant contribution from the circum- 



nuclear starburst, and the H2 emission we observe can be 
directly associated to this component. 

NGC 5728 also provides additional support to the idea 
of disconnection between the H2 emission and the NLR. 
This AGN is a classic example of a Seyfert 2 galaxy with 
a biconical emission line region separated by a dark band; 
the nucleus lies behind, at the point defined by the apexes 
of the two ionization cones (Wilson et al. I1993|l . The ion- 
ization cones are essentially collinear with the position 
angle of the brighter SE cone being 118°. In addition to 
the cones, the nuclear region of NGC 5728 shows other 
interesting features. In red and green optical continuum 
images presented by Wilson et al. 119931 , there are two 
miniature bars delineated by four main knots of emission, 
making the central region of NGC 5728 very complex with 
several decoupled dynamical components. In its K-band 
spectrum, NGC 5728 shows a remarkable II2 emission; it 
is the only object in which the II2 emission is stronger than 
Br7. Our spectrum covers the inner 320 pc with the slit 
oriented along the larger bar (PA=36°), almost perpen- 
dicular to the ionization cones and along the dark band. 
This suggest that the H2 emission is associated with the 
dark band, which serves as a natural reservoir of molecular 
gas. 

We also note that Mrk 1239 and H 1143-192 are the 
only Seyfert 1 galaxies of our sample with no II2 lines ob- 
served in their spectra. As [Fe 11] emission and other NLR 
features are clearly visible in their spectra (see Fig. ^, 
these objects provide additional support for a separate 
spatial origin between the II2 and NLR lines. 
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Table 3. Values of FWHM, in km s ^, for optical lines of 
the galaxy sample taken from the literature 



Galaxy 


[Oiii] A5007 


[Sii] A6717 


Mrkl210 


494^ 


276^ 


Mrk 1239 


630^ 


360"'^ 


NGC3227 


460^ 


427^ 


NGC4051 


162^ 


185^5 


NGC4151 


245^ 


280^ 


Mrk 766 


180^ 


1595 


NGC4748 


365' 


450'^,3 


NGC 5548 


411« 


300*^ 


Mrk 478 


920' 




PG 1448+273 


315' 




Mrk 493 


450' 




IH 1934-063 


562^ 


446-^ >3 


Mrk 896 


315' 




Ark 564 


265'' 


260* 



" FWHM of [N II] A6584. 

References — 1. Veron-Cetty et al. 2001; 2. Schulz & Henkel 
2003; 3. Rodriguez- Ardila et al. 2000; 4. Erkens et al. 1997; 5. 
Veilleux, 1991; 6. Moore et al. 1996. 



Finally from Table|21 the molecular and forbidden lines 
are all significantly narrower than the BLR lines. The last 
column of Table El shows the FWHM of Oi 1.128Aim, a 
pure BLR feature emitted in the outer portions of that 
region. The large difference in line widths between H2 and 
O I makes implausible the emission of H2 close to the BLR, 
i.e. from the torus of gas and dust around the central 
engine. 

3.1. NIR kinematics vs optical kinematics 

It is interesting to see how the widths of NIR lines mea- 
sured in §|21for the galaxy sample compares to those found 
in the optical region. This will allow us to get an insight 
of how the NLR gas is distributed and most importantly, 
to check if large amounts of high- velocity gas are obscured 
by dust extinction at shorter wavelengths. For this analy- 
sis we have chosen the lines [O ill] A5007 and [S 11] A6717. 
The former can be considered as the optical counterpart 
of [Siii] A9531 while the latter is representative of low- 
ionization gas as [Fe 11] is. Values of FWHM for these two 
lines, extracted from the literature, are listed in Table |31 
Note that for most cases, the spectral resolution is higher 
than that of NIR lines. 

Overall, Table|31supports the conclusions drawn in the 
previous section. Moreover, no significant broadening of 
NIR lines relative to the optical ones is detected. The ex- 
ceptions are Mrk 1210 and Mrk 1239. These two sources 
are known to have strong polarization due to dust scat- 
tering (Tran 1995; Goodrich 1989), supporting the view 
that high velocity gas are obscured by dust extinction. 
Note that for Mrk 766 and NGC 4051, the optical lines 
are half as broad as the NIR ones, but in both cases, the 
NLR lines are spectroscopically unresolved. 



It is worth commenting the results found for 
NGC 3227, NGC 4151, NGC 5548 and IH 1934-063. In 
these objects, the width of [Oili] is compatible to that 
found for [S III] , which argues against the presence of high- 
velocity gas in the visible region, not detected due to 
extinction. However, [Fell] is considerably broader than 
[S 11] (or [N 11] where the measurements of the sulfur line 
were not available). This gives additional support to the 
view that the [Fe 11] lines are at least partially formed in a 
distinct region from that originating other low-ionization 
species. Interestingly, the former three galaxies display 
radio-jets, whose interaction with the NLR has been pro- 
posed as the source of [Fen] lines. 

4. Excitation mechanisms of the NIR H2 lines 

The H2 molecule can be excited via three distinct mecha- 
nisms: (1) UV fluorescence, where photons with A < 912 
A are absorbed by the H2 molecule and then re-emitted 
resulting in the population of the various vibration- 
rotational levels of the ground electronic state (e.g. see 
Black & van Dishoek I1987|l : (2) shocks, where high- 
velocity gas motions heat, chemically alter, and accelerate 
the ambient gas resulting in excitation of the H2 molecule 
(e.g. see HoUenbach & McKee HM^ : and (3) X-ray illu- 
mination, where hard X-ray photons penetrate deep into 
molecular clouds, heating large amounts of molecular gas 
resulting in H2 emission (e.g. see Maloney et al. I1996|l . 
Each of these three mechanisms produces a different H2 
spectrum; therefore, the relative intensities between emis- 
sion lines of H2 can help to discriminate the dominant 
excitation process. 

Fluxes of the detected H2 and permitted Hi lines 
are listed in Table 0] They were measured by fitting a 
Gaussian to the observed profile and then integrating the 
flux under the curve. The LINER software (Pogge & Owen 
I1993|l was used for this purpose. For a large majority of 
our sample (16 out of 22 galaxies), these measurements 
are the first ones reported in the literature. H2(1,0)S(1) 
was detected at a 3a level in 20 out of 22 galaxies, the ex- 
ceptions being Mrk 1239 and H 1143-182 (see Figs.[T]toEl). 
It indicates that molecular gas within the inner 500 pc is 
common in AGN and not restricted to SB- AGN compos- 
ite objects. In addition, we clearly detected H2(1,0)S(2) 
for the large majority of AGN. We also detected strong 
H2(1,0)S(3) emission, sometimes larger than H2(l,0) S(l), 
but its flux must be used with caution because it may 
be severely affected by telluric absorption. Moreover, it is 
strongly blended with BrS to the blue and [Sivi] 1.963/im 
to the red. Fluxes for H2(1,0)S(0) and H2 (2,1)S(1) could 
also be measured at the 3cr level in 9 objects; fluxes were 
derived for the remaining objects using the 3 a error as an 
upper limit. In some Seyfert 1, it was possible to deblend 
the contribution of the NLR from the observed Br7 and 
Pa/3 lines. Their fluxes appears in Cols. 2 and 3, respec- 
tively, of Table El 

In order to see how the fluxes listed in Table 01 com- 
pare with previous measurements reported in the liter- 
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Table 4. Observed fluxes of H2, [Fell] and Hi lines in units of 10"^'^ erg cni~^ s^^. 





[Fen] 


Pa/Jt' 


[Foil] 


H2 


H2 


H2 


Br^t 


H2 


H2 


1.257Mm/ 


Galaxy 


1.257(/m 


1.282 Aim 


1.644 /jm 


1.957^111 


2.0332/jm 


2.1213Mm 


2.165 Mm 


2.247Aim 


2.223 Mm 


1.644Mm 


Mrkl210 


15.80±0.70 


31.45±0.88 


17.10±0.64 


3.01±0.30 


0.80±0.17 


2.00±0.12 


6.01±0.37 


0.44i0.09 


0.90i0.20 


0.92i0.05 


Mrkl24 


3.16±0.62 


14.6" 


2.82±0.47 


0.91±0.16 


0.36±0.05 


0.90±0.12 


2.49±0.55 


<0.18 


0.28i0.08 


1.12i0.29 


Mrkl239 


7.90±0.55 


198.6±2.60 


9.72±1.55 


<2.09 


<0.93 


<0.58 


32. Oil. 50 


<0.96 


<0.82 


0.81i0.14 


NGC3227 


41.1±1.18 


185.3±5.70 


41.0±3.60 


20.8±1.2 


7.97±0.92 


17.7±1.00 


30. Oil. 56 


2.32i0.21 


3.67i0.66 


1.00i0.09 


H 1143-192 


2.54±0.45 


182.0±2.77 


2.45±0.53 


<0.30 


<0.23 


<0.12 


26.7il.00 


0.24i0.08 


<0.14 


1.04i0.28 


NGC3310 


10.70±0.87 


31.74±1.03 


11.12±1.57 


<1.36 


1.24±0.25 


1.54±0.32 


10.79i0.33 


<0.32 


<0.67 


0.91i0.15 


NGC4051 


5.26±0.81 


88.4±1.80 


6.42±0.97 


7.51±0.62 


2.82±0.70 


5.81±0.37 


13.li0.84 


<0.56 


1.52i0.40 


0.82i0.18 


NGC4151 


58.8±2.94 


825.6±11.8 


56.2±1.96 


21.2±2.0 


6.84±0.42 


14.7±0.55 


153.8i7.7 


<0.62 


5.04i0.62 


1.04i0.06 


Mrk766 


7.69±0.70 


145. Oil. 89 


8.20±0.40 


4.39±0.39 


1.73±0.22 


2.36±0.24 


27.4i0.82 


0.30i0.06 


0.81i0.15 


0.94i0.10 


NGC4748 


7.39±0.37 


73.57±2.84 


7.85±0.61 


1.90±0.42 


<0.32 


1.59±0.20 


ll.32i0.67 


0.40i0.13 


<0.35 


0.94i0.09 


Mrk279 


7.22±0.46 


55.2±2.85 


5.36±0.61 


0.67" 


1.04±0.16 


1.89±0.23 


11.9i2.4 


0.56i0.15 


0.62i0.15 


1.35i0.18 


NGC5548 


1.71±0.26 


70.85±4.45 


1.30±a.l4 


<1.35 


<0.16 


0.80±0.11 


16.27i2.0 


<0.14 


<0.11 


1.31i0.24 


Mrk478 


<3.06 


82.70" 




3.73±0.40 


<0.90 


1.26±0.15 


6.90i0.48 




<0.35 




NGC 5728 


5.08±0.50 


7.37±1.16 


4.96±0.61 


6.04±0.56 


2.40±0.20 


6.28±0.12 


2.11i0.19 


0.89i0.20 


1.68i0.20 


1.02i0.16 


PG 1448-1-273 








0.93±0.07 


0.37±0.08 


0.45±0.08 


3.81i0.17 


<0.20 


<0.20 




Mrk291 


0.91±0.24 


4.27±0.44 


0.80±0.09 


0.30±0.12 


0.22±0.07 


0.47±0.10 


1.68i0.19 


<0.16 


<0.10 


1.14i0.33 


Mrk493 


1.96±0.34 


3.03±0.59 


1.4a±0.27 




0.32±0.07 


1.04±0.20 


3.45i0.25 


<0.26 


<0.21 


1.40i0.36 


PG 1612-1-261 


3.56±0.25 


50.92±0.98 




<1.02 


<0.07 


0.76±0.17 










Mrk504 


<0.72 


21.8±1.22 


<0.20 


<0.32 


0.15±0.04 


0.35±0.06 


3.48i0.45 


<0.13 


0.20i0.06 




IH 1934-063 


7.24±0.24 


159. Oil. 51 


7.46±1.05 


2.38±0.65 


1.80" 


1.23±0.20 


18.6i0.72 


0.31i0.10 


<0.10 


0.97i0.14 


Mrk896 


<0.40 


21.5±0.69 


<0.53 


0.41" 


0.42±0.13 


0.41±0.10 


3.21i0.59 


0.20i0.06 


<0.19 




Ark 564 


3.87±0.44 


90.9±1.50 


3.90±0.65 


0.80±0.22 


0.53±0.18 


1.24±0.27 


9.74i0.52 


<0.26 


<0.22 


0.99i0.20 



(a) Affected by telluric absorption. 

(b) Total flux of the line. 



Table 5. 

[Fen]. 



Fluxes of the narrow component of the permitted H i lines and line ratios measured between H i, H2 and 



Galaxy 



Br7 



Pa/3 



H2 2.12^m/Br7 [Fen] 1.257 fim/P&/3 



Mrk 1210 
Mrk 1239 
NGC 3227 
NGC 3310 
NGC 4151 
Mrk 766 
NGC 4748 
NGC 5548 
NGC 5728 
PG 1612+261 



6.01±0.37 
8.40±1.14 
6.85±2.00 
10.79±0.33 
23.26±1.49 
7.40±0.53 
1.39±0.26 
1.07±0.16 
2.11±0.19 
0.75±0.15 



31.45±0.88 
77.16±1.65 
23.34±1.78 
31.74±1.03 
102.4±2.52 
27.15±0.47 
8.77±0.50 
3.76±0.45 
7.37±1.16 
3.23±0.18 



0.33±0.03 

<0.07 
2.58±0.76 
0.14±0.03 
0.63±0.05 
0.32±0.04 
1.14±0.26 
0.74±0.15 
2.97±0.27 
1.01±0.30 



0.50±0.03 
O.lOiO.Ol 
1.76±0.14 
0.33±0.03 
0.57±0.03 
0.28±0.03 
0.84±0.06 
0.45±0.09 
0.69±0.13 
l.lOiO.lO 



ature, we have listed in Table El fluxes for H2 2.121/xm, 
[Fen] 1.64^m, Br7, [Fen] 1.257^m and Pa/3 published for 
a subsample of our objects along with the aperture size 
used in the observations. If more than one measurement 
is reported for the same object, they are listed in different 
lines. This is to avoid mixing data taken under different 
observing conditions and setup. Clearly, the effect of the 
extended emission is particularly important for NGC 4151 
and NGC 5728. For the remaining galaxies, most of the 
NIR emission seems to be concentrated in the circumnu- 
clear region. It important to keep in mind that the slit 
orientation may vary for each entry of Table El However, 
due to the relatively small apertures, its effects on the 
reported fluxes may not be significant. 

In the H band, we detected H2 (1-0)S(7) 1.747/im in 
NGC 5728 and NGC 3227; it was the only H2 line detected 
in that spectral region. This situation is in contrast with 
what is observed in planetary nebulae (Rudy et al. I2()()l|l 
and reflection nebulae (Martini et al. 11997(1 for example, 
where bright H2 lines from upper vibrational levels (v > 2) 
are usually present. In these sources UV fluorescence dom- 
inates the production of H2. In fact, the theoretical predic- 
tions of Black & van Dishoeck i(1987|l based on this mech- 



anism as the main source of energy input show that H2 
lines such as (6-4)Q(l) 1.601/im and (5-3)Q(l) 1.493/im 
may have fluxes comparable to that of (1,0)S(1). In our 
sample, the lack of molecular hydrogen from rotational 
levels with z^ > 1 in the H band suggests that UV fluo- 
rescence alone may have little importance in the produc- 
tion of the observed H2 spectrum. It is also possible that 
these features are diluted by the strong continuum from 
the AGN but this last hypothesis seems unlikely because 
[Fen] l.GAfim is clearly visible in all objects and it often 
has an intensity comparable to that of H2 (1,0)S(1). 

Studies aimed at investigating the origin of H2 lines 
carried out on AGN samples show controversial results. 
Veilleux et al. 1)1997(1 claim that shocks associated with 
nuclear outflows are a likely source of H2, while Quillen 
et al. p999.f) relate the H2 emission to local star forma- 
tion. Our data can be used to put stronger constraints 
on this issue by comparing the observed H2 fluxes with 
existing models that discriminate between these two pos- 
sible scenarios. Mouri 11994) proposed that the line ratio 
H2 (1,0)S(2) 2.247Aim/(l,0)S(l) 2. 2 12^m separates shock- 
dominated gas (with l2.247/im/l2.2i2/im=0. 1-0.2) from flu- 
orescent regions (with l2.247A«n/l2.2i2Mm ^0.55). However, 
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Table 6. Fluxes for the most important NIR lines taken from the literature. In units of 10^^^ erg cm^^ s^^ 



Galaxy 


H2 2.121^m 


[Fell] 1.64^m 


Br7 


[Fell] 1.257^m 


Pa/3 


Mrk 1210 


4.09' 




7.11' 




77.1' 


37.4' 


NGC3227 


29^ 

14±0.7'^ 
11.6±0.08* 


48±15^ 

32±1.7'' 

20.8±0.28* 


16.6±2.6^ 

16±1.1*^ 

22.1±0.18* 








NGC4051 


8.71 
4.2±0.5^ 


2.0±0.5'' 


13±1.9^ 








NGC4151 


2.0^ 








36.7±9'' 


836.2* 




50±5^ 


180±10'^ 


160±10^ 




230±10^ 


2300±50^ 


Mrk 766 


2.4I 












Mrk 478 






15.8±5^ 








NGC 5728 


2.6±0.4^ 














7.01^ 




2.22^^ 




29.9^ 


21.8^ 



References - 1. Ruiz 1997 (1".2 slit width); 2. Rudy et al. 2001. (3" slit width). 3. Schinnerer et al. 2001. (3". 6 circular aperture). 
4. Knop et al. 1996 (0". 75x1". 33 aperture). 5. Thompson 1995 (2"xl0" aperture). 6. Sossa-Brito et al. 2001 (2" aperture). 7. 
Veilleux et al. 1997 (3"x3" NGC 5728, 6"x3" Mrk 1210 in J; 3"x3" in K). 8. Reunanen et al. 2002 (l"xl".5). 



in dense gas (n > 10"* cm^^^) the coUisional de-excitation 
of the H2 molecule must be taken into account; the spec- 
trum is modified and approaches the thermal spectrum 
seen in shocks. 

Figure Ell shows (2,1)S(1) 2.247Aim/(l,0)S(l) 2.212/im 
plotted against (1,0)S(2) 2.033Aim/(l,0)S(0) 2.223/im. 
The models of Mouri (.1994|) are also shown. These lines 
ratios are reddening insensitive due to the very close wave- 
lengths. A large scatter is observed in both line ratios 
for the galaxy sample. NGC 3227, NGC 4051, NGC 4151 
and Mrk 766 (labelled 4, 7, 8 and 9, respectively) fall near 
the thermal curve model and indicate excitation temper- 
atures for the thermal component of between 1500 K and 
2500 K. For these galaxies, the H2 emission can be con- 
sidered purely thermal. For the remaining objects, a mix- 
ture of thermal and non-thermal processes are probably 
at work, with the excitation temperature of the thermal 
component higher than 1000 K. Interestingly, no AGN is 
located in the region occupied by the pure non-thermal 
UV excitation models, except possibly Mrk 896. However, 
its error bar makes any statement on this object highly un- 
certain. As expected, the starburst galaxy NGC 3310 falls 
far from the locus of points for the non-thermal models. 
For comparison, data for NGC 253 (Harrison et al. "1998"), 
a well-known starburst galaxy, are also plotted. NGC 253 
however, has a thermal excitation temperature below 1000 
K, similar only to Mrk 1210 and Mrk 504. Although it is 
not possible to state at this far what thermal mechanism 
excites the H2 lines, the results shown in Figure |S^ tell us 
that the thermal component leads over the non-thermal 
one. 

In order to quantify the contribution of the thermal ex- 
citation to the H2 lines, we have also plotted in Figure [51i 
the predicted line ratios from a mixture of thermal (model 
S2) and low-density fluorescence models (model 14) taken 
from Black & van Dishoek (1987). These two models were 
empirically mixed in varying proportions and the resultant 



emission line fluxes were then derived. The calculated line 
ratios are represented by the dashed line. It connects the 
models (represented by the triangles) where the percent- 
age of the thermal component decreases in steps of 10%, 
starting from a model where 90% is thermal and 10% UV 
fluorescence (the first triangle from left to right) up to a 
contribution of 20% thermal and 80% non-thermal. It can 
be seen that the observations that are scattered across 
the plot may, in a first approach, be explained by a suit- 
able contribution of thermal and non-thermal processes. 
Mrk 504 is probably the AGN with the largest contribu- 
tion of the later mechanism. In decreasing order of im- 
portance, we find that UV fluorescence responds for up 
to 50% of the observed H2 spectrum in Mrk 279, 30% in 
Mrk 1210, 25% in Mrk 124, and 15% in NGC 5728. 



Figure EJd shows H2 (2-l)S(l) 2.247Aim/(l,0)S(l) 
2.212/xm versus H2 (1-0)S(3) 1.963Aim/(l-0)S(l) 2.212/im; 
this is an additional diagnostic diagram proposed by 
Mouri H1994fl . equivalent to that of Figure |SJi. However, 
it has the drawback of including H2 (1-0)S(3) 1.957 ^m, 
which may be strongly affected by telluric absorption 
bands of H2O and CO2 if the redshift of the galaxy is 
between 2;=0. 017-0. 0398 or it may be strongly blended 
with [Sivi] 1.963/im. In addition, if the broad component 
of Br(5 is very conspicuous, it may affect the determina- 
tion of the continumn level around H2 (1-0)S(3) 1.957 /im. 
Overall, Figure \Sjp reproduces most of the conclusions al- 
ready drawn: NGC 4051 and NGC 3227 fall on or near the 
pure thermal excitation curve, and most objects show a 
thermal excitation temperature in the range 1000-3000 K 
indicating a dominance of thermal processes over the non- 
thermal ones. The departure of Mrk 766 and NGC 4151 
from the thermal curve can be explained by the reasons 
stated above. 



Alberto Rodriguez- Ardila et al.: H2 and [Fen] in AGN 



11 



CO 

CV2 



a. 

CO 
CO 

o 

C>j 



cvj 



OJ 



03 



2 - 



"T 1 T" 



n 1 T" 



"T 1 T" 




(a) 



-^^ 



15 



_L 



_L 



0.2 0.4 

Hg 2.247ij.m/2A2lfj.m 



0.6 



2.5 - 



T 



T 




(b) 



15 




_L 



_] I L. 



0.2 0.4 

Hg 2.247/im/2.12lMm 



0.6 



Fig. 5. Plot of (a) (2-l)S(l) 2.247Aim/(l-0)S(l) 2.121/xm vs (1-0)S(2) 2.033Aim/(l-0)S(0) 2.223^m and (b) (2-l)S(l) 
2.247Aini/(l-0)S(l) 2.121/xm vs (1-0)S(3) 1.957^ni/(l-0)S(l) 2.121/im . Curves represents thermal emission at 500- 
3000 K. Crosses are the non-thermal UV excitation models of Black & Van Dishoek H1987|l . Open circles are thermal 
UV excitation models of Sternberg & Dalgarno (,1989^1 . The numbers identify each object according to the notation 
given in column 1 of Tabled The black square in (a) represents the datum for NGC 253, a typical starburst galaxy, 

H1998|l . The full triangles, joined by a dashed-line, represents the predicted line ratios from 



taken from Harrison et al. 

a mixture of thermal and low-density fluorescence models of Black & Van Dishoek H1987fl . The percentage of the 
thermal component decreases in steps of 10% (and consequently the non-thermal component increases by the same 
proportion), starting from a model where 90% is thermal and 10% UV fluorescence (the first triangle from left to 
right) up to a contribution of 20% thermal and 80% non-thermal (last triangle of the sequence, to the right). 



4.1. Rotational and vibrational temperatures 

The main result found in the previous section can be con- 
firmed by deriving rotational and vibrational tempera- 
tures for the molecular H2 gas; thermal excitation must 
give similar rotational and vibrational temperatures, as 
would be expected for a gas in LTE. Fluorescent excita- 
tion, on the other hand, is characterized by a high vibra- 
tional temperature and a low rotational temperature; non- 
local UV photons, not characteristic of the local kinetic 



temperature, overpopulate the highest energy levels com- 
pared to that expected for a Maxwell-Boltzmann popula- 
tion. As a consequence, in photodissociation regions where 
UV fluorescence dominates, the measured vibrational tem- 
perature is ~5000K whereas the rotational temperature 
of the gas is ~10 per cent of that value (Sternberg & 
Dalgarno [TM^ . 

The values of T^ib and Trot found for the galaxy sam- 
ple are listed in Table [7| They were calculated using the 
fluxes of the observed H2 lines listed in Table ^ and the 
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expressions Tvib — 5600/ln(1.355 x /i-os(i)/-^2-is(i)) and 
Trot — -1113/ln(0.323x/i_os(2)/-^i-os(o)) from Reunanen 
et al. |2002). For Mrkl239, Mrk478, H 1143-192 and 
PG 1612-1-261, neither of the two temperatures could be 
derived either because only upper limits to the required 
fluxes are available or because the lines were not observed. 

Table strengths the hypothesis of the H2 lines be- 
ing purely thermal in NGC3227, NGC4051, NGC4151 
and Mrk766 because Tvib '^^^ Tlot. In Mrkl210, Mrk 504, 
NGC5728 and Mrk 896, the two temperatures differs ap- 
preciably, indicating the presence of non-thermal excita- 
tion, as was previously found from the sequence of mixed 
models of FigurelSJi. The values listed in TableQthen con- 
firms that for most objects, a mixture of thermal and non- 
thermal processes contributes to the observed H2 spec- 
trum. Except for Mrk 896 and Mrk 504, it can be stated 
that thermal processes are the leading excitation mech- 
anism, with a contribution that varies from 60% up to 
100%. 

The values for Tvib are also in very good agreement 
to those derived by Reunanen et al. l|2()()2|l for a small 
sample of AGN in which they found Tvib between 1800 K 
and 2700 K. For NGC3227, the only object common to 
both samples, a good agreement is found for Tvib with our 
results (2397±256K vs 1950 K of Reunanen et al. lM)2| . 

If the main excitation mechanism for H2 in AGN is 
thermal, it is interesting to examine the possible emis- 
sion mechanisms as either shocks (produced by supernovae 
winds or by a radio jet), or UV or X-ray heating are po- 
tential candidates for the II2 emission. 

Recently, Rodriguez- Ardila & Viegas H2()()3|l found ev- 
idence of starburst activity in Mrk 766, NGC4051, and 
NGC 3227 by observing the 3.3/xm polycyclic aromatic hy- 
drocarbon (PAH) emission in these objects. The beam size 
covered by their spectra is similar to the one used here. 
Galaxies with high levels of star formation activity fea- 
ture strong PAH emission and prominent H2 lines. The 
relationship between these two characteristics was initially 
established by Mouri et al. 1)1990(1 and later confirmed by 
Mizutani et al. H1994(l . who explained it in terms of the 
physical processes in photodissociation regions associated 
with relatively massive young stars. Our data provides fur- 
ther support to the hypothesis of the starburst as the ori- 
gin of the observed molecular hydrogen emission in those 
three objects. In this scenario, shocks driven by the super- 
nova remnants appears as the most likely source for the 
H2 excitation. In an AGN with a starburst component, 
part or all of the observed molecular emission is expected 
to arise from this process. This seems to be the case of 
Mrk 766, NGC 4051, and NGC 3227. The predicted H2 2- 
1S'(0)/1-0S'(1) ratio for the pure shock model (S2) of Black 
& van Dishoek (1987) at T=2000 K is 0.082, in marginal 
agreement to the measured ratios of 0.13±0.03, 0.096, and 
0.13±0.01, respectively. Their S2 model also predicts a H2 
l-05'(0)/l-0S'(l) ratio of 0.21, weU in accord to the values 
of 0.34±0.07, 0.26±0.07, and 0.21±0.04, respectively. 

Though supernova-driven shocks offer a plausible sce- 
nario to explain the pure thermal H2 spectrum observed 



in those three objects, we cannot ignore that AGN are 
powerful sources of X-ray emission. It is then expected 
that X-ray heating may have a non-negligible contribu- 
tion to the observed H2 flux. Moreover, supernova rem- 
nants themselves emits high-energy radiation capable of 
heating the gas. 

For NGC 4151, the non-detection of the 3.3 /im feature 
within the inner 3". 8x3. "8 region (Imanishi et al. 19^?H)l 
makes unlikely the hypothesis of a starburst as the main 
contributor of the H2 lines. Further efforts to detect sig- 
natures of starburst activity in that Seyfert have failed 
(e.g. Sturm et al. I1999|l . The analysis of the molecular 
line profiles carried out in § |3 showed that it is unlikely 
that H2 is produced by shocks from the radio jet due to 
the small value of FWHM measured. Thus X-ray heating 
is the most plausible mechanism for H2 excitation in this 
object; specific models to test this possibility are discussed 
in§0 

Finally, we comment on Mrk 1210. For this galaxy, we 
derived a Trot ~ 900 K. We already know that UV fluo- 
rescence has at most, a 30% contribution. It means that 
the low Trot can be attributed to UV heating and this pro- 
cess can be the main source of H2 excitation. UV heating 
models (Black & van Dishoek Il987|l . where FUV photons 
emitted by OB stars (photodissociation regions) heat the 
gas, predict excitation temperatures of <1000K for the H2 
thermal component. Such low values are typically found 
in starburst galaxies (Mouri 1994|, suggesting that the H2 
emission observed in Mrk 1210 may primarily come from a 
circumnuclear starburst. In fact, the position of Mrk 1210 
in Figure Ot is almost coincident with that of NGC 253, 
a well-known starburst galaxy. Additional evidence giv- 
ing support to this hypothesis was found in Section 
where the kinematics of the H2 gas pointed out that it 
was unrelated to the AGN. Further arguments supporting 
this view come from Storchi-Bergmann et al. 1)1998(1 and 
Schulz & Henkcl ( 2003| who reported the detection of a 
Wolf-Rayet (WR) feature around A4686 A in the inner 2" 
of Mrk 1210. This aperture coincides with the beam size of 
our observation. Moreover, WR-features are taken as un- 
ambiguous signs for the presence of a starburst (Schulz & 
Henkel U()()3(l and due to the close relationship between H2 
and starburst activity (Mouri 02201 Mizutani et al. 11994(1 
it is natural to attribute the H2 emission to this last com- 
ponent. Observations at higher spatial resolution can test 
this scenario. 

In summary, the observational evidence presented up 
to now shows that the molecular hydrogen lines in AGN 
are mostly thermal, with a contribution of non-thermal 
processes typically varying from 15% to 30%. There is not 
unique process associated with the thermal component. 
H2 gas heated by UV photons emitted from photodisso- 
ciation regions. X-rays from the AGN, and shocks from 
supernovae remnants arise as the most likely sources of 
H2 excitation. Our results agree with those of Rigopoulou 
et al. 1)200 21. found using mid-infrared H2 lines from ISO 
observations for a sample of Seyfert 2 galaxies. 
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Tvib 


J rot 




Tvih 


J- rot 


Galaxy 


(K) 


(K) 


Galaxy 


(K) 


(K) 


Mrkl210 


3100±700 


900±300 


NGC 5548 


<2700 




Mrk 124 


<2900 


2500±800 


NGC 5728 


2500±600 


1400±400 


NGC3227 


2400±300 


3100±700 


PC 1448+273 


<5000 


<2200 


NGC3310 


<3000 


<2200 


Mrk 291 


<4100 


<3300 


NGC4051 


<2100 


2200±800 


Mrk 493 


<3300 


<1600 


NGC4151 


<1600 


1300±200 


Mrk 504 


<4300 


800±300 


Mrk 766 


2400±500 


3000±700 


IH 1934-063 


3300±1200 




NGC 4748 


3300±1200 




Mrk 896 


5500±2100 


<3300 


Mrk 279 


3700±1100 


1800±500 


Ark 564 


<3000 


<4400 



5. The origin of the [Feii] Lines 

The origin of [Fen] in active galaxies is highly contro- 
versial. It can be produced by shock excitation from the 
radio jets and/or photoionization from the central source. 
In addition, fast shocks associated with supernova rem- 
nants can be an additional source of [Fen], as is fairly 
well established in starburst galaxies (Simpson et al. 119961 
Moorwood & Oliva 1988 ). As found by Forbes & Ward 
P993|l . the strength of the [Fen] emission in the central 
regions of AGN is tightly correlated with the 6-cm radio 
emission, and that the same correlation is obeyed for both 
starburst and Seyfert galaxies. This was taken as evidence 
that the most likely mechanism for the production of [Fe n] 
emission in Seyfert galaxies is shock excitation from the 
radio jets and/or SNRs. 

From Table 01 and Figures ^ to |2| it can be seen that 
[Fen] 1.257/im and [Fen] 1.64 /im were detected in all ob- 
jects of our sample but PG 1448+273. In the latter AGN, 
both lines are located in regions of very poor atmospheric 
transmission, hampering any effort to detect them. 

The flux ratios [Fen] 1.644/im/Br7, or equivalently, 
[Fen] 1.257/xm/Pa/3^ can be used to study the origin of 
the [Fe n] emission because there is an increasing progres- 
sion in the ratio from pure photoionization (as in Hn 
regions) to pure shock excitation (i.e. supernovae rem- 
nants) , with starburst and active galaxies located at inter- 
mediate values (Mouri et al. 119931 Goodrich et al. 119941 
Alonso-Herrero et al. ll997|l . In starburst galaxies the ratio 
[Fen] 1.644/im/Br7 is observed to be in the range 0.5—2 
(Simpson et al. [I996), in reasonable agreement with the 
starburst models of Colina (1993) which predicts ratios 
between 0.1 and 1.4. In these objects, the [Fen] emission 
originates in supernova-driven shocks (Mouri et al. UOOO ) . 
AGN usually show larger values of [Fen]/Br7 than star- 



^ The flux ratios [Fell] 1.644/im/Br7 and [Fell] 
1.257/ini/Pa/3 are equivalent because the two [Fell] lines 
share the same upper level. They are related by the equa- 
tion [Fell] 1.257Mm/Pa^ = 0.231*[Fell] 1.644^m/Br7. The 
advantage of using [Feii] 1.257/im/Pa/3 instead of [Feii] 
1.644/im/Br7 is that the former is reddening insensitive due 
to the proximity in wavelength of both lines, serving as a more 
reliable diagnostic. 



burst galaxies, probably reflecting the effect of the Seyfert 
nuclei (Kawara et al. ll988|l . In fact. X-ray emission, which 
is dominant in Seyferts, can penetrate deeply into atomic 
gas and create extended partly ionized regions where [Fe n] 
can be formed. Models presented by Alonso-Herrero et al. 
II1997|I show that X-rays are able to explain [Fen]/Br7 
ratios up to ^^20. 

Col. 5 of TableElhsts the flux ratio [Fe n] 1.257 fim/Paf] 
as derived for our galaxy sample. Only the objects in 
which it was possible to carry out a clean deblending of 
the narrow component of the permitted lines are included 
in this analysis. Overall, our line ratios agree with the 
values found for other samples of Seyfert galaxies, 0.5 < 
[Fe n] 1.257//m/Pa/? < 4.6 (see Fig. 3 of Alonso-Herrero et 
al. I1997|l . Interestingly, Mrk 766, a NLSl with a circum- 
nuclear starburst, has a [Fen] 1.257/im/Pa/3 ratio simi- 
lar to that observed in starburst galaxies, indicating that 
most of [Fe n] we detect is of stellar origin. On the other 
hand, NGC 3227, a classical Seyfert 1 also harboring a 
circumnuclear starburst, has a high [Fen] 1.257/im/Pa/3 
ratio, pointing towards an AGN origin. NGC 4151 and 
NGC 5548 are in the borderline of the starburst- AGN divi- 
sion but any contribution from star forming regions would 
be negligible because of the lack of circumnuclear star- 
burst in these objects. 

Although the [Fen] 1.257/im/Pa/3 is helpful to discrim- 
inate between a stellar and non-stellar origin for [Fe n] , it 
can say little about the excitation mechanisms that can 
lead to the [Fen] lines in the latter case. Heating either 
by X-rays or shocks created by mass outflows from the 
nuclei (i.e., jets and/or winds interactions with ambient 
clouds) are plausible scenarios. As proposed by Mouri 
et al. H2000|l . these two processes can be discriminated 
by the electron temperature of the [Fen] region: 8000 K 
in heating by X-rays, and 6000 K in heating by shocks. 
Determining the [Fen] temperature, however, is a diffi- 
cult task, mostly because some of the diagnostic lines are 
intrinsically weak and heavily blended with broad adja- 
cent features. The alternative is to examine if X-rays can 
explain the observed emission. This issue is addressed in 
detail in Section l5^ 
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Table 8. Line ratios H2 2.12r7zum/Br7 and [Fell] 1.257/im/Pa/3, from AGN taken from the literature. 



Galaxy 


Ha 2.121^m/Br7 


[Fell] 1.257Aim/Pa/3 


Reference 


NGC 1386 


1.14±0.10 


1.59±0.14 


1 


NGC4945 


3.13±0.06 


0.86±0.03 


1 


NGC 5128 


2.04±0.10 


3.42±0.16 


1 


NGC 1097 


>4.40 


>0.43 


1 


Mrk 1066 


0.89±0.06 


0.75±0.04 


2 


NGC 2110 


3.17±0.63 


8.11±0.78 


2 


NGC 4388 


0.88±0.08 


0.40±0.03 


2 


Mrk 3 


0.31±0.04 


1.24±0.04 


2 



References: (1) Reunanen et al. f^KMl : (2) Knop et. al. lOIMl 



5.1. Are [Fen] and H2 related? 

The observational evidence presented throughout this 
work indicates that [Fell] and H2 are common features 
in the nuclear spectra of AGN, although probably orig- 
inated from different parcels of gas. However, this does 
not exclude the possibility that both set of lines originate 
from a single, dominant mechanism. Moreover, if dust is 
intermixed with the line emitting clouds, and there are 
strong velocity gradients, [Fe 11] and H2 lines can appear 
to have different velocity fields (Knop et al. I2001|l even 
though they are produced in adjacent regions. 

In order to examine a possible relationship between 
[Fell] and H2, we have plotted in Fig.ElH2 2.121/xm/Br7 
versus [Fe 11] 1.257^m/Pa/3 for our sample along with data 
taken from the literature for Seyfert 2 galaxies observed 
with a similar beam size and resolution (see Tabled . Note 
that the line ratios are completely reddening free due to 
their close proximity in wavelength for the lines involved. 

Interesting results are seen in Fig. El For most ob- 
jects with both Il2/Br7 and [Feii]/Pa/3 <2, the two ratios 
shows a possible, positive correlation. However, for values 
of either ratio larger than 2, the correlation breaks down 
and no trend is observed. In addition, no Seyfert 1 galaxy 
displays extreme values of any of the line ratios plotted 
in Figure El Both H2/Br7 and [Fell] /Pa/3 >2 occurs pre- 
dominantly in Seyfert 2 galaxies. 

The correlation observed in Fig. El is consistent with 
the one found by Larkin et al. (,1998 , ) in a study of LINERs 
and other emission line objects, including a few Seyferts. 
They report a strong linear correlation in the log- log space 
between [Feii]/Pa/? and H2/Br7, with LINERs display- 
ing the largest ratios and starbust galaxies the smallest 
ones. Larkin et al. I|1998(l suggested that [Feii]/Pa/3 ^1 
and H2/Br7 ~3 marked the end of Seyfert- like nuclei and 
the beginning of LINER-like objects. The larger sample 
of Seyfert galaxies employed here allow us to propose that 
the borderline between Seyfert and LINER activity occurs 
at [Feii]/Pa/3 --2 and H2/Br7 -2. In this sense. Fig. El 
would serve as a diagnostic diagram in the NIR to sepa- 
rate galaxies according to their level of nuclear activity: 
starburts occupying the bottom left corner of the plot, 
with [Feii]/Pa/3 and H2/Br7 <0.3, Seyferts occupying the 
locus of points with values of either ratio between 0.4-2, 



and LINERS showing ratios larger than 2. We note that, 
as in our plot, Larkin et al. (1998) found that some of 
the Seyfert 2 galaxies lie significantly below or above the 
correlation. 

Assuming that [Fen] /Pa/3 and H2/Br7 are in fact cor- 
related, interpreting the origin of that correlation is non- 
trivial. This is mainly because Fig.Eldoes not discriminate 
between the different processes that may give rise to ei- 
ther [Fell] or H2 in each object. The most straightforward 
interpretation is that both set of lines are excited by the 
same mechanism. It is yet possible that more than one 
mechanism drives the production of [Fe il] and H2 within 
the same object, but environmental effects dilutes the ef- 
fects of one mechanism over the other (i.e., the presence 
of a circumnuclear starburst in an AGN). 

Recalling that the objects plotted in Fig. El are all 
AGNs, it is natural to consider that the dominant excita- 
tion mechanism discussed above is related to the central 
engine. In particular. X-rays from the nucleus can heat 
the NLR gas through the large partially ionized zone it 
creates, driving the [Fe 11] and H2 emission. Since 98% of 
the iron is tied up in dust grains, this process must free up 
iron through dust destruction and yet not destroy the H2 
molecules. Blietz et al. ()1994fl and Knop et al H2UU1|I show 
that hard X-rays can do it for a sample of Seyfert 2 galax- 
ies. Additional evidence of the role of X-rays can be drawn 
from the good correlation between H2 2.121/xm/Br7 and 
[O i] 6300A/Ha and between [Fe 11] 1.257/iim/Pa/3 and [O l] 
6300 A/Ha reported by Larkin et al. ifTO^ . They indeed 
argue that hard X-ray heating from a power-law source 
is a plausible mechanism to explain LINERS with low 
values of the [Feii]/Pa/3 and H2 2.121/im/Br7 ratio, to 
which they suggest, are low-luminosity Seyfert galaxies. 
This same picture can be extended to the Seyfert galax- 
ies of our sample to support the slim correlation reported 
here. 

5.2. X-rays as a source of [Feu] and H2 lines 

In order to test if the observed [Fe 11] and H2 emission can 
be attributed to X-ray heating, we have calculated the 
emergent H2 2.121/im and [Fell] 1.644/im fluxes using the 
models of X-ray excitation of Maloney et al. fl996). These 
models predict the resultant infrared emission (lines and 
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0.5 -i—i- 



Mrk 3 







A 

NGC 1097 
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log(H,/Br7) 



0.5 



Fig. 6. The ratios H2 2.121fim/BTj and [Fell] 
1.257/im/Pa/3. Only the Seyfert Is and NLSl in which 
a rehable deblending of the narrow component of the 
Hi Hnes was carried out were plotted. Filled circles 
represents our data and open symbols data taken from 
the literature. Triangles are measurements from Knop 
et al. H2U01|I and squares from Reunanen et al. I|2U02|I : 
see text for further details. The number along each data 
point identifies our galaxies according to the notation 
given in column 1 of Table Q 

continuum) of a neutral gas cloud of electron density ric 
exposed to a source of hard X-rays of brightness L^ lo- 
cated at a distance d from the cloud. The main parameter 
controlling the physical state of the X-ray-illuminated gas 
is the effective ionization parameter, Ccff , which is the at- 
tenuated X-ray flux to density ratio, determined from the 
expression 



c 



,eff 



100 



L44 



n^dl^N^i 



(1) 



where L44 is the hard X-ray luminosity, L-^, in units of 
10** erg sec~^, n^ is the gas (total hydrogen) density in 
units of 10^ cm~^, dpc is the distance, in parsecs, to the 
AGN from the emitting gas, and N22 is the attenuating 
column density between the AGN and the emitting gas, in 
units of 10^^ cm^^. Once Ceff is computed, the predicted 
H2 and [Fell] emergent intensities can be obtained with 
the help of Figs. 6a and 6b of Maloney et al. H1996|l for a 
gas density of 10^ cm""^ and 10'^ cm~'^, respectively. 

The input parameters employed in our modeling are 
listed in Cols. 2-4 of Table El Two values of the hydro- 
gen column density, Nh, are listed and were used in the 



modeling: one determined directly from X-rays (Col. 3) 
and one from 21 cm radio observation (Col. 4). If this last 
value was not available, we used instead the Njj derived 
from E(B-V) according to the expression Hhi = 5.2 x 
10^1 E(B-V) cm-2 (ShuU & van Steenberg [MSI , which 
assumes Galactic dust/gas ratio. The distance d between 
the X-ray source and H i cloud was assumed to be equal to 
the radius of the integrated region (see Col 10 of Table^. 
This value represents the maximum distance at which the 
line could be emitted. 

As a first test, we determined the predicted H2 
2.121/im intensity for a gas density nc=10^ cm~'^. The 
results appears in Cols. 6 and 7 of Table O for Nh de- 
rived from X-rays or either radio observations or E(B-V), 
respectively. Compared to the observations, H2 intensi- 
ties calculated using Nh obtained from X-ray data offer a 
poorer fit. This can be understood if we take into account 
that the hydrogen column density derived from hard X- 
rays maps material obscuring the innermost part of active 
galaxies, while the H2 and [Fen] emission is expected to 
arise farther out in the NLR. For this reason, we opted to 
use the Nh values determined either from 21 cm measure- 
ments or E(B-V). With these constraints, we calculated 
the emergent [Fe 11] and H2 intensities and the results are 
listed in Table [TOI 

The results presented in Table nTH suggest that X-ray 
heating is responsible for only a fraction of the [Fe 11] and 
H2 emission, except in NGC 5548 where it can fully drive 
it. In four out of seven objects, the predicted [Fell] flux is 
neghgible while in NGC 4151 and NGC 3227, X-rays may 
contribute for up to 15% of the measured [Fe 11] . Regarding 
II2, X-rays contribute to the observed emission but only 
in NGC 4151 NGC 5548, and NGC 3227 can it entirely ex- 
plain the excitation of the molecular gas. Note that for the 
Seyfert 2 galaxies Mrk 1210 and NGC 5728, an additional 
source of H2 would be necessary because X-rays contribute 
with only a small fraction of the measured flux. 

The above discrepancy can be alleviated if the emitting 
clouds were effectively located closer to the X-ray source 
than the conservative distance we adopted in the models 
(equals to the beam size of our observations). In Mrk 1210, 
for instance, neutral clouds with ric=10^ cm"'^ at ~80 pc 
can generate enough II2 and [Fe 11] to account for the ob- 
servations. A similar result is obtained for NGC 5728. This 
possibility would obviate the need of additional excitation 
mechanisms and readily explain the correlation observed 
in Fig. El In order to confirm or discard this scenario, data 
taken at a superior resolution in order to map the inner 
molecular gas and constrain its density and distance from 
the central source are necessary. For the time being, we 
are only in a position to state that overall, the role of 
X-ray heating in exciting H2 and [Fen] lines in AGN is 
confirmed but the strength at which it operates strongly 
varies from object to object. 

Morphological effects, more visible in Seyfert 2 than 
in Seyfert 1 galaxies, can be behind the strong departure 
from the correlation of Fig. presented by some objects 
of the sample. In effect, most of the galaxies with large 
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Table 9. Predicted H2 2.12/j,m fluxes from tfie X-ray models of Maloney et al. H1996|l for a density of 10^ cm~'^. 





log Li 


iVHx 


A'h.nlr 




Predicted H2 2.12/im 


Galaxy 


(erg s 1) 


(10^1 cm-2) 


(10^1 cm-2) 


p2 


A^Hx 


-^H NLR 


(1) 


(2) 


(3) 


(4) 


(5) 


(6) 


(7) 


Mrkl210 


42.36 


0.83 


0.286'' 


2.00±0.12 


0.32 


0.58 


NGC3227 


42.00 


0.50 


0.036 


17.7±1.00 


1.29 


0.68 


NGC4051 


41.20 


6.30 


0.156^ 


5.81±0.37 


0.16 


1.51 


NGC4151 


43.00 


2.26 


0.380 


14.7±0.55 


1.45 


24.6 


Mrk766 


43.00 


1.28 


0.260^ 


2.36±0.24 


0.67 


1.68 


NGC5548 


43.76 


4.57 


0.021 


0.80±0.11 


0.74 


283.0 


NGC 5728 


41.34 




0.391^ 


6.28±0.12 




0.19 



Notes. Col(l) — Object name, Col. (2) — Logarithmic 2-10 keV luminosity, Col. (3) — hydrogen column density determined from 

X-ray observations. Col. (4) — hydrogen column density for the NLR determined either from radio 21 cm observations or the 

E(B-V). Col. (5) — Observed H2 2.12/im flux. Cols. (6-7) — predicted H2 2.12/im flux using column density from X-ray and radio, 

respectively. 

1 ASCA 2-10 keV luminosities were taken from George et al. HM^ for NGC 3227, NGC 4051, NGC 4151, Mrk 776 and NGC 5548 

and from Awaki et al. JM^Hll for Mrk 1210. 

^ In units of 10~^^ erg s~^ cm~^ . 

^ Determined from the relation Nhi=5.2x10^^E(B-V). See text for further details. 

Table 10. X-ray model parameters and predicted line fluxes for [Fell] 1.64/im and H2 2.12/im. 





Tic 

(cm-3) 


(10^1 cm-2) 


[Fell] 1.64 /im" 


H2 2.12 


/im^ 


Galaxy 


Obs. 


Predicted 


Obs. 


Predicted 


Mrk 1210 


10^ 


0.286 


17.1±0.64 


0.0 


2.00±0.12 


0.58 


NGC 3227 


10^ 


0.036 


41.0±3.60 


7.08 


17.7±1.00 


33.9 


NGC 4051 


10^ 


0.156 


6.42±0.97 


0.0 


5.81±0.37 


1.51 


NGC 4151 


10^ 


0.380 


56.2±1.96 


10.7 


14.7±0.55 


24.6 


Mrk 766 


lO'^ 


0.260 


8.20±0.40 


0.0 


2.36±0.24 


1.68 


NGC 5548 


10^ 


0.021 


1.30±0.14 


15.8 


0.80±0.11 


10.96 


NGC 5728 


10^ 


0.391 


4.96±0.91 


0.0 


6.28±0.12 


0.19 



In units of 10 



erg s 



H2 2.121/im/Br7 or [Fell] 1.257/im/Pa/3 ratios share sim- 
ilar morphological aspects. For example, the presence of 
ionization cones and, most importantly, dust lanes, favour 
the production of H2 over [Fe 11] . This could be the case in 
NGC 4945 and NGC 5728, where such features are clearly 
detected (Reunanen et al. I^OH^ Wilson et al. ITO^ . In 
these objects, H2 2.121/im/Br7 >2. Jets and outflows, on 
the other hand, tend to favour [Fell] emission over H2. 
This is the case in the galaxies NGC 2110 (which shows 
one of the [Feii]/Pa/9 ratios ever measured in an AGN, 
Knop et al. '"2001), NGC 5128 (Centaurus A) , and Mrk3. 
In this case, [Fell] 1.257/im/Pa/3 >2. Whether these fea- 
tures are in fact related to the enhancement of any of the 
two emissions cannot be deduced from our data. Higher 
spatial resolution is needed to address this question. 

Finally, it is interesting to note that X-rays models 
can explain [Fcii]/Br7 '^20, as found for Alonso-Herrcro 
et al. H1997|l , but it would need quite a high X-ray flux or 
clouds located too close to the central source. None of the 
objects studied (except probably NGC 5548), has the right 
combination of X-ray emission, hydrogen column density, 
and distance from the X-ray source to the emitting clouds, 
to produce ratios as high as 20 (using Maloney 's models). 



Since these variables were constrained from observations, 
the results reflect the conditions for the galaxies of our 
sample. 



5.3. The [Few] lines as useful indicators of reddening 
for the NLR 

The simultaneous measurements of [Fell] 1.257/im and 
1.644/im allow us to test, for the first time on such a big 
sample of AGN, how theory compares to observation. The 
ratio [Fe 11] 1.257/im/ 1.644/im is fixed by atomic physics to 
be 1.34 (Bautista & Pradhan ll998|l since both lines orig- 
inate from the same upper level (i.e. '^Dt/2-^D'j/2I'^D7/2- 
^-^9/2)- It means that this ratio is not only a good check 
on theoretical predictions but also a reliable indicator of 
the extinction towards the NLR because it is indepen- 
dent of temperature and density. Last column of Table 0] 
lists the [Fell] 1.257/im/ 1.644/im ratio measured in each 
galaxy. All but Mrk 279, NGC 5548 and Mrk 493, tend to 
have ratios ^30% lower than the theoretical prediction. 
In fact, the mean [Fell] 1.257/im/ 1.644/im of the sample 
equals 0.98. The uncertainties in the observed line ratios 
are 10%-20%. 
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The above result can be interpreted in two ways. Either 
the the bulk of [Fell] arises in a dusty region, or the A- 
values for the transitions involved need a fine tuning. In 
the first case, it implies that the extinction affecting the 
gas emitting the [Fe il] lines is large and surprisingly ho- 
mogeneous. Assuming an intrinsic 1.257/xm/1.644/im ratio 
of 1.34, the observed ratios agree with Av ~3.7 (adopt- 
ing the Cardelli et al. 1 19891 extinction law and foreground 
screen). This last value would be even higher in Mrk 1239 
and NGC 4051 but overall, rather similar in most AGN. 
Note that the extinction towards the NLR determined by 
other authors for some of the objects of our sample is sig- 
nificantly lower: Av ~0.6 for NGC 3227 (Crenshaw et al. 
EOnTJ), ~0.43 for Ark 564 (Crenshaw et al. lMI^ and -1.27 
for Mrk 766 (Gonzalez Delgado & Perez [TMSJI . 

In the second case, our data can be useful to set im- 
portant constraints to a more precise determination of the 
4-values for [Fell] forbidden lines. Bautista & Pradhan 
l|1998jl had already noted about the difficulty associated to 
the determination of accurate forbidden transition prob- 
abilities, especially for Fell, owing to the large numbers 
of algebraic terms involved. Unfortunately, there are very 
few works published in the literature with a systematic de- 
termination of the ratio /(1.257/im/1.644/^m) on a similar 
sample. However, our results tend to agree to that of Sugai 
et al. (dnnnil, who reports a ratio of 0.80 for ARP 299. 
The extinction Av, determined by Heisler & De Roberties 
H1999|l from /(1.257/im/1.644/j,m) on a small sample of Hii 
and Seyfert 2 galaxies, also tends to be significantly higher 
that that derived from optical data. 

In short, our observations imply that the extinction 
measured from the [Fell] lines is large compared to that 
determined from optical indicators. This can be due to a 
genuine effect and it is somehow expected because the NIR 
probes a larger optical depth than the optical. An alter- 
native is to consider that the intrinsic 1.257/xm/1.644/im 
ratio is ~30% smaller than the actual theoretical value. 
This is certainly an issue that deserves further attention 
and need to be explored using a larger sample of objects. 

6. Final Remarks 

In order to discuss the origin of the molecular hydrogen 
and [Fe 11] emission lines, near infrared observations have 
been carried out for a sample of 13 NLSl, 6 Seyfert 1, 2 
Seyfert 2 and 1 starburst galaxies. In the following, our 
main conclusions are briefly reviewed. 

The kinematical link between the molecular hydrogen 
and the NLR gas was analysed on the basis of the line 
widths. For all galaxies in the sample, H2 2.121/im was 
in the limit of the spectral resolution (360 km s~^) and 
narrower than the NLR lines for the cases in which the 
latter lines were resolved. [Fell] was resolved for 50% of 
the sample. In NGC 3227, NGC 4151 and PC 1612+261 it 
was broader than [Svi] 1.963^m. It is then unlikely that 
II2 and [Fe 11] come from the same volume of gas. 

The detection of H2 (1,0)S(1) 2.121/imin 20 out of the 
22 objects of our sample reveals, for the first time, strong 



evidence that molecular gas within 500 pc from the cen- 
tre is common in AGN and not restricted to SB/AGN 
composite galaxies. H2 lines from vibrational transitions 
were detected but are not as bright as in planetary nebu- 
lae, providing the first indication that UV fiuorescence is 
not the main mechanism that excites them. This result is 
strengthened by the analysis of the line ratios that discrim- 
inates between shock-dominated and fiuorescence produc- 
tion. For 4 objects (NGC 3227, NGC 4051, NGC 4151 and 
Mrk 766) the excitation mechanism is clearly thermal. For 
the remaining AGN, a mixing with a non-thermal process 
cannot be discarded, although the results point to a dom- 
inant thermal mechanism. This result is confirmed by the 
similarity between the vibrational and rotational temper- 
atures of the II2. However, a unique source of the thermal 
component could not be established. The relative contri- 
bution to this component coming from shocks originating 
in SN explosions or in radio jets, and from X-rays heating 
may vary from one object to the other. 

The origin of the [Fell] lines is also controversial. [Fell] 
1.257/xm was detected in all galaxies of our sample, sug- 
gesting a hnk to the AGN itself. The [Fen] 1.257/xm/Pa/3 
ratio is used to study the origin of the [Fe 11] emission. It 
reveals that stellar and non-stellar processes contributes 
to the [Fell] spectrum with a predominance of the lat- 
ter in most objects. The [Fen] 1.257/j,m/Pa/3 ratio values 
are within the interval usually observed in other AGN. 
Mrk 766 is particular in this respect, however. Being a 
NLSl, the observed [Fen] is compatible with a starburst 
origin. 

A correlation between H2 2.121/im/Br7 and [Fen] 
1.257/xm/Pa/3 was found but it breaks down for the galax- 
ies where either one of the ratios is larger than 2, pre- 
dominantly Seyfert 2s. Starburst galaxies are preferen- 
tially located in the region with [Fe n]/Pa/3 and H2/Br7 < 
0.4 while LINERS are characterized by [Fen] /Pa/3 and 
H2/Br7 > 2. We confirm the usefulness of this plot as a 
diagnostic tool at separating emitting line objects by ac- 
tivity level and set constraints to the locus of points shown 
by AGN. 

Line intensities predicted by X-ray heating models are 
able to explain the H2 and part of the [Fen] emission 
in Seyfert 1 but fail at reproducing these characteristics 
in Seyfert 2. The results are improved if the emitting 
clouds are located closer to the central source, —80 pc 
or less. This would explain the correlation between H2 
2.121/im/Br7 and [Fen] 1.257^m/Pa/3 discussed above. 
However, the presence of significant starburst activity in 
some of the objects makes it plausible that more than 
one exciting mechanism is at work. In Seyfert 2, X-rays 
appears to have a lesser importance and the H2 emission 
is probably enhanced by contributions from circumnuclear 
starburst and structures such as ionization cones and dust 
lanes. Nonetheless, the role of these structures is not con- 
firmed by the present data. 

Finally, our data allowed us to study, for the first time, 
the [Fe n] 1.257/im/1.644/im ratio in Seyfert 1 galaxies. We 
found that most objects have values 30% lower than the 
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intrinsic value predicted by theory. It implies either than 
the extinction towards the [Fell] emitting clouds is very 
similar in most objects and significantly larger than that 
derived from optical observations, or there is an overesti- 
mation of the ^-values for the Fell transitions. We dis- 
carded a systematic error on the measurements because 
the uncertainties associated to the observed ratio is be- 
tween 10%-20%. 
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